o 
o 

(N 
< 



00 

o 



Oh: 



X 



ZZ^ and Z77 couplings at linear e+e collider energies with the effects of Z 

polarization and initial state radiation 

S. AtafI and i. gahiii] 

Department of Physies, Faculty of Sciences, Ankara University, 06100 Tandogan, Ankara, Turkey 

The constraints on the neutral gauge boson couphngs, ZZy and ^77, are investigated at linear 
e''"e~ collider energies through the Z'y production with longitudinal and transverse polarization 
states of the final Z boson. Because of high energy of linear electron- positron beams, initial state 
radiation (ISR) considerably changes the production cross section. We obtain an increase in the 
cross section by a factor of 2-3 due to ISR for transverse polarization and by a factor of 10-100 for 
longitudinal polarization states depending on the energy. 

We find the 95% C.L. limits on the CP conserving form factors /if, /if, and hj with integrated 
luminosity 500/&~^ and =0.5, 1, 1.5 TeV energies. It is shown that the longitudinal polarization 
of the Z boson, together with ISR, can improve sensitivities by factors 2-3. 

PACS numbers: 12.15.Ji, 12.15.-y, 12.60.Cn, 14.80.Cp 



I. INTRODUCTION 



The self interactions of the gauge bosons are consequences of the non-abelian structure of the electroweak sector of 
the Standard Model (SM). The study of the trilinear gauge boson couplings leads to important tests of electroweak 
J> [ interactions. Neutral gauge boson couplings ZZ^ and Z77 are not generated at tree level by the SM. Higher order loop 
CO I level corrections are higly below the current experimental sensitivity 0. Nevertheless, the new physics with energy 
scale above present experimental threshold might provide tree level neutral triple-gauge boson couplings. Deviation 
of the couplings from the expected values would indicate the existence of new physics beyond the SM. Therefore 
precision measurements of triple vector boson vertices will be the crucial tests of the structure of the SM. 
, For the process e+e^ — > Zj it is convenient to study the anomalous neutral triple gauge couplings ZjZ* and ^'77* 
' which obey Lorentz and gauge invariance. Within the formalism of Ref. P, Q there are eight anomalous coupling 
parameters hf , HJ {i = 1, .., 4) which are all zero in the standard model. Here we are interested in CP-even couplings 



that are proportional to and {V = Z,j). The photon and the Z boson in the final state are considered as 



3 

1 on-shell particles while the third boson at the vertex, the s-channel internal propagator, is of-shell. Due to partial 
wave unitarity constraints at high energies, an energy dependent form factor ansatz can be cosidered: 

(1 + S/A2)3 



hi i^) = TTT^Tij^ ; *-i,3 (1) 



hns)^7—\^ ; * = 2,4 (2) 



C3 ■ . . . 

■ ■ ' In this work we assume that new physics scale A is above the collision energy y/s and we neglect the energy dependence 

of the form factors in the energy region we are intereted in. 

CP conserving anomalous Z(pi)^{p2)Z[p^) vertex function can be written following the low energy parametrization 
of the residual effect from the effective lagrangian P, Q : 



^9eTzyZ[PiiP'2iP3) = We 2 



"-3 e P2p + ^P3 e I 



(3) 



where mz and ge are the Z-boson mass and charge of the proton. The Z77 vertex function can be obtained with the 
replacements: 



2 9 2 

P3-P1 _^ Pi_ 

2 2 ' 



/jf^/i7, * = 3,4 (4) 
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The overall factor in the ^77* vertex function originates from electromagnetic gauge invariance. Due to Bose 
statistics the ^77 vertex vanishes identically if both photons are on shell (Yang's theorem) Q . 

Previous limits on the Z7Z* and ^77* anomalous couplings have been provided by the Tevatron |/if | < 0.36, 
\hf\ < 0.05, \m < 0.37 and \hl\ < 0.05 H and the combination of four LEP experiments ALEPH, DELPHI, L3, 
OPAL -0.20 < /if < 0.07, -0.05 < /if < 0.12, -0.049 < hj < 0.008 and -0.002 < hj < 0.034 ^] at 95% C.L. . 
Based on the analysis of ZZ production at the upgraded Fermilab Tevatron and the CERN Large Hadron Collider 
(LHC) achievable hmits on the ZZ7 couplings ( ff, fj, f^, /g) have been discussed 0,0]. 

Research and development on linear e+e^ colliders at SLAC, DESY and KEK have been progressing and the physics 
potential of these future machines is under intensive study. In this paper, the process e^e~ — > Zj with the final state 
i~^£~'y is investigated to search for ZjZ* and .^77* couplings. Because of the high energy of the incomig beams, 
initial state electromagnetic radiative corrections (ISR) are taken into account using the structure function method. 
Another important point is the polarization of the final state Z boson. In order to determine the Z polarization, we 
show that the angular distribution of the Z decay products, has a clear correlation with the helicity states of the Z 
boson. 



II. CROSS SECTIONS AND ANGULAR CORRELATIONS FOR FINAL STATE FERMIONS 

In this section, we present the cross section calculation via helicity amplitudes for the complete process e+e^ — > 
Z"f — > £^£^^ and describe angular distributions of final state fermions to see the correlations with the polarization 
states of Z boson. Let us start with the differential cross section 



where pi , p2 are the momenta of incoming leptons, p^ , p4 are the momenta of outgoing fermions and k2 is the momentum 
of outgoing photon. |Mp is the square of the full amplitude which is averaged over initial spins and summed over 
final spins. The helicity dependent full amplitude can be expressed as follows 



{2n)H\pi + P2 - fci - fc2)(2^)4(54(fci - p3 - p4) (6) 

where fci is the internal momentum of the Z boson propagator, cti, (T2 are the incoming lepton helicities, (73, (74, are 
the outgoing fermion and photon helicities. For the summation over the intermediate Z boson polarization we take 
the helicity basis Xz — +,—,0. Here Dz{ki) is the Breit-Wigner propagator factor for the Z boson 

= kl-^4 + ^Mz^z 

Here AIa{(7i, a2', Xz , X^) is the helicity amplitudes for e+e^ — > Z7 with on shell Z boson which are provided in the 
appendix. M;,(A2; cs, (74) is the decay amplitude of the Z boson which will be given later. Using the narrow width 
approximation differential cross section becomes 

da{e+e- ^ Z7 - i+rj) = ^p^^l E Ma{Xz)AMXz)\' 

d^ki d^k2 ,^ 44 



d^ps d^pi 
(27r)32£;3 {2'KY2Ei 



{2nY5''{pi+P2-ki-k2) 



(2^)V(fci-p3-P4) (8) 



Here Ma{Xz) and M}j{Xz) indicate average over initial lepton spins, summmation over final state fermion spins and 
photon. In the rest frame of the Z boson where the decay amplitudes are most simply expressed, we get 
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Az 

-(27r) (5 (pi +P2 - Ki - ^2) 



d cos e*d(j)* (9) 

where 6** and are the polar and azimuthal angles of the final state leptons in the Z rest frame with respect to the 
Z boson direction in the ££j rest frame. After integration over azimuthal angles (j)* interference terms will vanis as 
below 

da{e+e- ^ ^ i+r^) = l^^^-2__[|Af^(+)M,(+)|2 + \Ma{-)Mt{-)\' + |M.(0)M,(0)|2] 

d^ki d^k2 44 



This can be written in a more compact form 



(27r)*r(pi +P2 - fci - fc2)dcosr (10) 



d<j{e+e- ^ Z7 ^ i+r-i) = . ^ . [\daii+)\Mb{+)\^ + dai{-)\Mi,{-)\^ + dai(0)|Afb(0)|2] dcos^ (11) 

[62TrMzi z) 

If one performs further integration over polar angle cos 9* it turns out to be a well known result 

da{e+e- ^ Z-f ^ e+Tj) = [dcri(+) + dcri(-) + dai{0)]BR{Z -> £+r) 

= dai{e+e- ^ Zj)BR{Z ^ e+r) (12) 

Explicit forms of the decay amplitudes |Af6(+)p, |Mb(— )p and |Mb(0)p in the Z rest frame are given by 

=^[ffi(l-cosr)2+g|,(l + cosr)2] (13) 

\MtH\' =^[5i(l + cosr)2+g2^(l-cosr)2] (14) 

|M,(0)|2 = MU9I+ 9%) sin^ 6* (15) 
\M,iTR)\^ = |M,(+)|2 + |Mb(-)|2 (16) 
|A4(L0)p =|Mb(0)p (17) 

By measuring the polar angle distributions of the Z decay products one can directly determine the differential cross 
sections for fixed Z helicities. In other words, Z helicity states are obtained from a fit to these distributions. Complete 
factors {1 /32Trr z)\Mi,{Xz)\'^ in front of d(Ti{Xz) in differential cross sections are plotted in FigQ] As can be seen from 
Figniongitudinal(LO) and transverse (TR) distributions are well separated from each other. This is why we consider 
only transverse and longitudinal polarizations for the Z boson in next sections. 



III. INITIAL STATE ELECTROMAGNETIC RADIATIVE CORRECTION 

Due to small mass of the electron, a significant role is played by the electromagnetic radiative corrections to the 
initial electron-positron state epecially at linear collider energies. In this work we use structure function formalism 
to describe the electromagnetic radiative corrections in e~^e~ colliders 't'I . The cross section can be written in the 
following form within this formalism 



a{s) ^ dxi dx2 Di{xuQ^) D2{x2,Q^y{s') 



(18) 
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FIG. 1: Polar angle distributions of Z decay product {l/'i2-KV z)\Mb{\z)\^ in the rest frame of Z boson for \z = —1,-1-1,0. 
Transverse polarization state is defined as the sum of -f-1 and -1 states. 



where (t'{s') is the cross section with reduced energy s' = xiX2S. Di{xi,Q'^) {D2{x2,Q'^)) stands for the electron 
(positron) structure function giving the probability of finding an electron (positron) within an electron (positron) 
with a longitudinal momentum fraction xi {x2)- Although several definitions of the structute functions are present 
we use the following ones which are used by HERWIG 



Di{x,s) ^ P{l^xf~^g{x,s) 



g{x,s) 



f3{l+x/2)x/2 



0' 



4 log X 



(1 - P'-) + y[(l + x){{l + xY + Slogo:} - 



(log^-1), y^m-xf-^] 



(19) 
(20) 
(21) 



To avoid divergency at the upper limit of the momentum fraction, x = 1 , the cross section can be transformed into 
different form 



a(,s) 







l-€ pl-e 

dxi I dx2 Di{xi,s) D2{x2,s)a'{s') 



l-e 



+e'^g2{l,s) dx2 Di{xi,s)a'{xis) + e'^gi{l,s) dx2 D2{x2, s)a' {x2s) 



+e^^gi{l,s)g2{l,sy{s) 



(22) 



where e can be taken as 10^^ — 10^^^. In this region of e the cross section changes by a factor of 0.985. If one takes 
smaller e values, higher machine precision gives softer e dependence. The following transformation gives relatively 
smooth integrand 



dxi / dx2 Di{xi,s) D2{x2,s)(t'{s') = / dE dF gi{xi, s) g2{x2, s)(t' {xiX2s) (23) 



where 
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XI = 1 - i-Py/^ , X2 = l- {-EY'P (24) 

Emin = —(1 — j^™™ ) , -Emaa; = (25) 

F ■ - (I 'I P - -f^ ('261 

man — V / ; max — c V ^/ 

- M r27i 

S 

Using above formalism we calculate the cross section with initial state radiative corrections for both Standard 
Model and anomalous coupling cases. In Fig|21the effects of ISR to the total cross section for the standard model 
process e+e" are shown as a function of energy y/s (before ISR) with transverse(TR) and longitudinal(LO) Z 

polarization. Here average over initial spins and sum over photon polarization are performed. The unpolarized cross 
section is almost the same as the TR polarization case since the magnitude of the cross secion is dominated by the 
cross section with TR polarization. As can be seen from the figures, ISR increases the cross section with a factor 2-3 
for the TR case and a factor of 10-100 for LO polarization case depending on the energy. Furhtermore, very small 
LO cross section becomes sizable due to ISR. The reason for this comes from the fact that the energy dependences of 
the LO and TR cross sections are different. In order to understand this feature let us write standard model squared 
amplitudes for transverse and longitudinal Z polarizations using the helicity amplitudes given in the appendix 

\M^{TR)\'^ 8[(C3^)2-f(C3«)2](-^-l)|^i±^ (28) 



sm 



\Ma{LO)\'^ 16[{Ctf + {C^rh ^ (29) 

where Ma{Xz) is defined in Eq.(8) and we neglect electron mass. As can be seen here the longitudinal part is 
independent of the polar angle whereas the transverse part strongly depends on it. Angular integration implies that 
contribution of the transverse polarization always gives larger cross section than the longitudinal one. In both cases, 
the largest cross section comes from the energy region where ^/s ~ mz- When we consider ISR, the major contribution 
from the integration over xi and X2 to the cross section is due to the lower limit of the s' = xiX2S 
limit the above amplitudes take the forms 



\Ma{TR)\^^ l&[{C!^f + [Cif][^-l)-^ (30) 

sm {iXsy 

4 

\AULOr^ lQ[[C^f + [Cif]^ (31) 

where we use s = m?^ -f As with As << m?^. 

The cross section without ISR leads to s' = s with xi = X2 = 1. For the collider energy y^s = 1 TeV where 
s >> m\ the amplitudes becomes 

|M„(Ti?)p^ S[{C^f + {Cif]{^-l){l + 2^) (32) 

sm 3 

|M,(iO)p~ 16[(C|')2 + (C3«)2]I^ (33) 

s 

Here the LO cross section continiously decreases as s increases. This is the expected result for longitudinal polarization. 
Because there is no LO polarization for masless vector bosons. For y/s = lOmz the effect of ISR on both amplitudes 
can be compared easily as below: 

\MaiTR)\jsR ^ 2m| ^^^^ 



\Ma{TR)\^ - (As) 
\Ma{LO)\jsH ^ 2m 



|M,(iO)|2 - (As) 



50--^ (35) 
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FIG. 2: Effect of ISR to the standard model cross section for e e "fZ with transverse and longitudinal polarization of Z 
boson. Two polarization states, TR and LO, get different contributions from ISR. 



FigOshows the energy dependence of the total cross sections with anomalous coupling parameter /if = 0.01 for TR 
and LO polarization states. The contributions of anomalous couplings become remarkably important after the center 
of mass energy 500-600 GeV. Effect of ISR is also shown for two cases. In the case of LO polarization, after the energy 
of 900 GeV the ISR gives negligible corrections. Similar behaviour appears for hj = 0.01 and hf — hj — 0.001 values. 
In Fig^ the energy dependences of the ISR corrected total cross sections with four anomalous coupling parameters 
are plotted for LO polarization and unpolarized cases. 

It is also important to see how the anomalous couplings change the shape of the angular distribution of the Z boson 
for the polarized and unpolarized cases. Angular distribution of the Z boson is shown in Fig. [Slfor LO polarization 
state with and wihout ISR correction. Since the TR polarization state of the Z boson is poorly sensitive to anomalous 
couplings, angular distributions with TR case are not plotted. In all figures, only one of the coupling parameters are 
kept different from zero. From all these figures we reach at the following remarkable results. The LO polarization 
states are always more sensitive to the anomalous couplings. Much larger deviations arise from hX for both TR 
and LO polarizations. Z'yy* couplings always provide the higher contribution to the cross section than the Z^Z* 
couplings. The shape of the curves differs for two kinds of polarizations of the Z boson. The ISR correction gives 
larger contribution to the SM cross section than the case with anomalous couplings. Therefore, the sensitivity of the 

e~ 7Z process to the anomalous couplings is expected to become poorer due to the ISR correction. Numerical 
results for all polarization configurations will be given in the next section. 



IV. LIMITS ON THE ANOMALOUS COUPLING PARAMETERS 



If the Z boson decays into a pair of charged leptons the signal for the final states can be ^11 where we consider 
t — e,ii. The potential background processes for £ = /i final state are the following: 

e+e^ Z{'j) jii : s channel Z or 7 exchange (final state bremstrahlung). 

e^e~ — > 77 j££ : t channel e exchange. 

For £ = e final state additional t channel background processes such as those arising from both Z(7) and e exchange; 
Z or 7 exchange (final state bremstrahlung) are present. Since we take into account only on shell Z bosons, we should 
impose a cut on the invariant mass of charged leptons Mu — Mz- This cut reduces the effect of background processes 
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FIG. 3: Energy dependence of the total cross section for e'^e~ jZ with anomalous ZZ'y vertex parameter h§ = 0.01. 
Effects of ISR and Z polarization are also shown. As in the previous figure, two polarization states of the Z boson get different 
contributions from ISR. 



foi £ — ^ drastically. The total cross section of background processes at least 100 times smaller than the process 
e'^e^ Z'j (t channel e exchange with on shell Z boson). In the case of ^ = e final state background cross 

sections are 10 times higher when compared to /i final states. Therefore the major potential background is due to 
£ = e final states. In the following sensitivity calculation, background contributions have a negligible effect. In order 
to obtain realistic limits on the h\ and h\ from the linear e^e'' collider the number of events have been calculated 
using N = Aa{e'^e~ ^Z)BR{Z £~^£~)Lint for integrated luminosity Lint = 500/6~^. Here lepton channel of 
the Z decay and overall acceptance A=0.85 has been taken into account. For total cross section | cos9\ = 0.99 has 
been used as angular region. The 95% confidence level (C.L.) limits have been estimated from total cross section 
using simple one parameter test for =0.5, 1, 1.5 TeV. At ^/s = 1.5 TeV the SM cross section without ISR 
correction gives the smallest number of events N=40 (using the above formula) for LO polarization state. If ISR 
correction is included number of events increases up to N=2800. For lower initial energy, -^i = 1 TeV, the smallest 
event number becomes 250 and 10000 without and with ISR correction. The limits which have been obtained are 
shown in Tables IHIIII for the deviation of the cross section from the standard model value without systematic error. 
It should be noted that better limits are obtained for the polarization configuration Xz =L0 which leads to the order 
of 0(10-3 _ 10-4) for /if and hj, 0(10-^ - lO'^) for /if and hj at ^ = 0.5 - 1.5 TeV when ISR corrections are 
taken into account. Without ISR correction, the LO polarization of the Z boson improves the limits by factors 3-7 
depending on the energy. The LO polarization with ISR correction improves the limits by factors 2-2.5 which should 
be taken as the realistic results. As can be seen from the tables, TR polarizations are not sensitive to anomalous 
couplings. But the advantage of TR polarization case is the absence of /if and hj couplings. This feature reduces the 
number of coupling parameters. In order to see the degree of energy dependence on the anomalous couplings let us 
take into account the increase in cm. energy from 0.5 TeV to 1.5 TeV for the LO polarization configuration. Then 
we get the improvements in sensitivity limits by a factor 20, 30 for and /ig, by a factor 150, 200 for hf and hJ. 



V. CONCLUSION 



Future linear e+e" colliders with = 0.5 — 1.5 TeV energy and integrated luminosity 500/6-"'^ probe the ZjZ* 
and Zjj* anomalous couplings with far better sensitivity than the present colliders Fermilab Tevatron and LEP2 
experiments. Measurement of final state Z boson polarization is important in two ways. First, LO polarization state 
is always far more sensitive to anomalous couplings. Second, TR polarization state contains only and hJ. Inital 
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FIG. 4: Energy dependence of the total cross section for e^e^ "fZ with LO polarization (left) and unpolarized (right) cases. 
The values of the anomalous ZZ'y and Z'yy vertex parameters hs = hj = 0.01 and hf = hj = 0.001 are chosen. All cross 
sections are ISR corrected. 
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FIG. 5: Angular distributions of the outgoing Z boson in the center of mass frame of outgoing Z and the photon. 6 is the angle 
between the outgoing Z boson and the incoming electron. Angular dependence of the four anomalous couplings are shown 
with (left) and without(right) ISR corrections. Only LO polarization states are taken into account in these figures since TR 
polarization states are poorly sensitive to anomalous couplings. 
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TABLE I: Sensitivity of the linear e^e collider to ZZ-y and ^77 couplings at 95% C.L. for — 0.5 TeV and Lint = 500 
/6"\ Only one of the couplings is assumed to deviate from the SM at a time. 
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|3 X IQ-^I 
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TABLE II: Sensitivity of the linear e+e" collider to ZZ-y and Z77 couplings at 95% C.L. for ^ = 1 TeV and Li„t = 500 /6~\ 
Only one of the couplings is assumed to deviate from the SM at a time. 
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|3 X 10-^1 
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state electromagnetic radiative correction improves the number of events especially for the Standard Model processes 
at the linear collider energies. Furthermore, LO polarization case gets larger contribution from the ISR than the case 
of TR polarization. 

Some limits on the above couplings via e~^e^ —^ Zj process at linear collider energies from similar works are the 
following for comparison: the sensitivity (one standard deviation) 2 x 10~^ , 4 x 10"'^, 4 x 10~^, 3 x 10~^ for hJ, 
/if, hJ, /if, respectively at 500 GeV energy and 100 luminosity 0; the sensitivity(95% C.L.) 0(10-2) for hJ, 
/if, 0(10-3) for hJ, /if at 500 GeV and 10 fb''^ Predictions from CERN LHC || give the limits 5.2 x lO^^, 

3.7 X 10-^ for /if and hf at 2a, 10 fb~^. Comparison of the results in this work at y/s = 1.5 TeV with those of LHC 
shows that our results are improved one order of magnitude. For more precise results, further analysis needs to be 
supplemented with a more detailed knowledge of the experimental conditions. 
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APPENDIX: HELICITY AMPLITUDES 



There are four Feynman diagrams for the process e+e" if one includes Z-yZ*, Z77* vertices. HeUcity 

amplitTidcs Mi and M2 arc responsible for the diagrams concerning Z^jZ* and .^77* interactions arising from s 
channel Z or 7 exchanges. My, and M^ are standard model contribution of the t and u channel of the process. The 
parameters of the helicity amplitudes M(a-e, cr^; Az, A-y) are helicities of incoming electron and positron, outgoing Z 
boson and photon. The values they take are given by : 

(Te:L,R , u'^:L,R , Az:+,-,0 , A^ : +, - (A.l) 

Here L and R stand for left and right. Helicity amplitudes we have obtained for each Feynman diagram in the cm. 
frame of e+e" can be written as follows: 



M{ae,a'^;Xz,X'r) = Mi(CTe, cr^; Az, A^) + M2{ae,(j'^;Xz,X'r) 

+M3((7e, a'^; Xz, A^) + M^^e, u'^; Xz, X^) (A.2) 



Mi{LR; 
Mi{RL; 
Mi{LR; 
Mi{RL; 
Ml (Li?; 
Mi{RL; 
Ml (Li?; 
Mi{RL; 



++) 
++) 
+-) 



-) 
-+) 
-+) 



Ml (Li?; 0-t-) 
Mi(i?L;0+) 
Mi(Li?;0-) 
Mi(i?L;0-) 



-C^hliMl 



s) sin 6 
s) smO 



= -C^hiiMl 
= 
= 

= Cfiif(M|-s) sin 6* 
= -Cf/if(M|-s)sin6' 
= 
= 



^/2Mz 



Vs(M| - s)(l +COS6') 



V2Mz 
V2Mz 



V2M, 



ys(M|-s)(l-cos6') 
^/s{Ml- s){l- cose) 

y/s{Ml- S){1 + cose) 



{Ml s)hf 



2M| 



s)hl 



2M| 



(M|. 



s)hl 



2M| 



+ 



(M| 



s)hl 



2M| 



(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.10) 

(A.11) 
(A.12) 
(A.13) 
(A.14) 



^ 2M| ' 1 2M| 



gL = ^{Cv + Ca) , 9R = ^{Cv - Ca) 

Cv = 2sm^ew-l , Ca = -\ 

gz = . — ^ , gl = 47ra 
sin 0w cos ffw 



M2{LR 


++) 


= -C2hl(Ml - s) sin 61 








M2{RL 


++) 


= M2{LR-++) 








M2{LR 


+-) 


= 








M2{RL 


+-) 


= 








M2{LR 


-+) 


= 








M2{RL 


-+) 


= 








M2{LR 




= C2hl{Ml~ s)siTie 








M2{RL 




= M2{LR; ) 








M2{LR; 0+) 


= -^^Vi(M|-s)(l + cos0) 


hj- 


{Ml-s)hl 

^ 2M| 


M2{RL;0+) 






2M| 


M2{LR- 0-) 


= M2{RL;0+) 








M2(-RL;0-) 


= M2{LR; 0+) 
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r,7?- -i—u] 


Cfs(l-cos6l) sing 
2Af| + (s - Af|)(l - cos e) 


(A.32) 


Mi( RTr -\—^-^ 


C§[2Ml - s(l - cosg)] sing 
2M2 + (s-Af|)(l-cosg) 






^71-8(1- cos g) sing 

2M| + (.s-Af|)(l-cosg) 


('A 34"! 




Cfs(l + cosg) sing 
2A^2_^(s-Af|)(l-cosg) 


(A.35) 


M-i(LR- — ) 


C^[2Ml- s{l-cose)]sme 
2Af| + (s-Af|)(l-cosg) 


(A.36) 




Cf s(l -cosg) sing 
2A//2 + (s-A//2)(l-cosg) 


(A 371 


MadR: — h) 

J. r-i. Q I / 1-'^ 1 / 


C|^s(l + cosg) sing 
2A^2 _^ _ Ml){l - cos g) 


(A.38) 


M^{RL--+) 


r';f .s(l - corigj siu^^ 
~ 2Af| + (s - Af2)(i _ cosg) 


(A.39) 


M3{LR;0+) 


V2Mz[2M^ + (s - A'/|)(l - cosg)] 


(A.40) 


M3{RL;0+) 


C3^y2i[3A^2 -s+(s + Ml) cos g] sin^ | 
Afz [2A^2 + (s _ M| ) ( 1 - cos g)] 


(A.41) 


MsiLR; 0-) 


C| v^[3i\/| - ^ + (/V + cos g] siu^ f 
Afz[2Af| + (s - A//2)(l - cosg)] 


(A.42) 


M3(i?L;0-) 


C^^s{s + Af|) sin^ g 
y2Afz[2Af| + (s - Ai"2 )(1 - cos g)] 


(A.43) 



where 

= Qe9e9L , Ci = QegeGR (A.44) 
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Mi{LR;++ 
Mi{RL-++ 
Mi{LR-+- 
Mi{RL;+- 

Mi{LR- 

Mi{RL- — 
Ma{LR--+ 
Mi{RL--+ 
Mi{LR; 0+ 
Mi{RL; 0+ 
Mi{LR; 0- 
Mi{RL; 0- 



Cf [s(l + cos6i) - 2M|]sin6' 



2M| + (s - 


M2)(1 +cos6») 


C^s{l 4 


- cos 9) sin 6* 


2M2 + (s - 


M2)(1+cos6») 


Cf s(l - 


- cos 6) sin ^ 


2M| + (s - 


M2)(l + cos6') 


Ci«s(l 


+ cos 9) sin 9 


2M2 + (s 


-M2)(1 + cos6') 




- cos 6*) sin 9 


2M| + (s - 


M|)(l + cos6») 


Cf [s(l + cos 61) - 2M|] sine* 


2M2 + (s - 


Af2)(i +cos6») 


Cfs(l 


+ cos 61) sin9 


2M2 + (s 


- Mj)(l + cos9) 




- cori 0) sill 6^ 



2M| + (s-M2)(1 + cos0) 
Cf y2i[s + (s + Ml) cos6' - 3M2] cos 
Mz[2M| + (s - M2)(1 + cos6»)] 
C^y^js + Ml) sin^ 61 



2 e 

2 



x/2Mz[2M| + (s - M|)(1 + cos^)] 
Cf^/^Cs + MDsin^g 
\/2Mz[2M| + (s - M|)(l + cos 6*)] 
Ci^v^[s + (s + M|) cos 9 - 3M|] cos 
Mz[2M2 + (s - M|)(l + cos 61)] 



(A.45) 
(A.46) 
(A.47) 
(A.48) 
(A.49) 
(A.50) 
(A.51) 
(A.52) 
(A.53) 
(A.54) 
(A.55) 
(A.56) 



where 



i — ^3 ^ 



(A.57) 



